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CN ' Abstract 

o. 

00 , We propose an inclusive analysis of a stransverse mass {rriT2) using a hemisphere 

method for supersymmetry studies at the LHC . The hemisphere method is an 
algorithm to group collinear and high pt particles and jets, assuming that there are 
two of such groups in a event. The rriTi is defined as a function of the unknown 
LSP mass, two hemisphere momenta, and missing transverse momentum. The 
kinematical end point of the niT2 distribution provides information on the squark 
and gluino masses. We perform a Monte Carlo simulation to study the inclusive 
m,T2 distribution at the LHC. We show that the end point of the inclusive mn 
distribution has a cusp- like structure around the true LSP mass. The knowledge of 
the expected kinematical behavior near the end point for true events is important 
to establish the end point of the inclusive distribution. We hnd that the inclusive 
analysis is useful to obtain the information on the heaviest of the squark/gluino. 



1 Introduction 



While the particle interactions at low energy are described correctly by the standard 
model (SM), the mechanism of electroweak symmetry breaking by Higgs boson suffers 
from the fine turning problem. In addition, the SM is not successful to describe the dark 
matter in our Universe. 

We expect to obtain information on the physics beyond the standard model at ATLAS 
and CMS experiments at the LHC, which is scheduled to start this year (2008). Among 
the various proposals, the phenomenology of the models with quark and gauge partners 
with multiplicatively conserved parity, such as supersymmetric models with conserved R 
parity, Little Higgs models with T parity and and universal extra dimension models, get 
much attention. In the supersymmetric models, quark and gluon partners (squark and 
gluino) are pair produced at the LHC, and subsequently decay into the SM particles and 
the lightest supersymmetric particles (LSP's). The signature at the LHC will be the high 
Pt jets and leptons with significant missing transverse momentum which arises from the 
LSP escaping from detection. By the end of the LHC experiment, the squark and gluino 
in minimal supergravity model will be searched up to ~ 2.5 TeV [TJ [2]. 

Interests in the new physics go beyond the discovery. Many studies have been carried 
out to find out the possible clues to study "the nature of the new physics", such as the 
masses, spins and interactions of the new particles. The progress has been made especially 
in the exclusive channels. The end points of the invariant mass distributions constrain the 
sparticle masses and for some cases nearly all sparticle masses can be measured. The end 
point study is extremely successful when the decay involves many leptons [U [2j [3J HJ [5J [6] . 

Other important variables are transverse masses. The peak of the effective mass 
distribution, which is the sum of the transverse momenta of jets, leptons and E™ ss , 
is related to the sum of the sparticle masses produced by pp collision. The quantity is 
inclusive and would be useful in the early stage of the LHC experiment. More sophisticated 
quantity is the m<r 2 variable [TJ [8]. This can be calculated from the two visible objects, 
and the missing momentum of the events and a test LSP mass. Recently, the quantity 
acquires much of attention because this function has a cusp at the correct LSP mass when 
the squark/gluino undergoes three body decays. Some exercises have been carried out for 
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some model points where only gluino - gluino production can be observed, or for several 
other decay patterns without specifying the selection processes [9j [TOj [TTJ [12]. 

In this paper, we propose an inclusive study of variables using a hemisphere 
method. Namely we group jets into two "visible objects" and calculate rriT2 variable 
based on them. The grouping algorithm is called a hemisphere method and discussed 
in earlier works [131 EE]- The motivation is to collect the cascade decay products from a 
squark or a gluino enough probability to see the rriT2 end point, and obtain the information 
on their masses without going into exclusive analysis. If it works, the LSP mass also can 
be obtained in the early stage of the experiment. 

We found the correspondence between the reconstructed rriT2 end point and the mass 
of squark or gluino is good. We recognized that m^2 is sensitive to max(m§, trig) as 
rriT2 is defined from the maximum of transverse mass of the visible objects for the test 
LSP momenta consistent with E™ 1SS . Not surprisingly, the probability to reconstruct the 
correct object is rather low and the end point is smeared, however, we find the event-wise 
response of the m,T2 to the test LSP mass mentioned in Refs. [91 [10J, [TH [12] is useful to 
ensure the correctness of the end point. 

In this paper, we especially compare the mixed modulus anomaly mediation (MMAM) 
model to the supergravity (SUGRA) model. The MMAM predicts a degenerate mass 
spectrum in some parameter region of the model, where the sparticles are heavy while 
available p^s of the daughter particles are small. This is the model where the kink at 
the LSP mass should appear in a rather high value. This point may be compared with 
model points in the mSUGRA (minimal Supergravity) model, where the gluino mass is 
lighter than the corresponding MMAM point, but the squark mass is much heavier than 
the gluino so that the total squark/gluino production cross sections are same. We find 
that the end point is useful to extract the squark mass for this case, therefore the 
MMAM and the SUGRA models can be distinguished. 

This paper is organized as follows. In Section [2] we review the rriT2 variable and the 
cusp structure appearing in the endpoint of the m<r 2 distribution as the function of the 
test LSP mass. In Section [3] we describe the inclusive and perform a Monte Carlo 
simulation to study the distribution. Section H] is devoted to the conclusion. 
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2 Transverse mass (m^) 



In hadron collisions, squarks and gluinos are produced in pair and these SUSY particles 
decay subsequently into the final states including jets, leptons, and two LSP's. The LSP 
is usually the lightest neutralino. With i?-parity conservation, the LSP is a neutral and 
stable particle and it escapes from detection. There are two LSP's in the final state 
and one cannot measure each LSP momentum experimentally while the total transverse 
momentum can be measured. The stransverse mass rriT2 is defined as follows: 



m T2 {m x ) = mm 



max^p^.p^ 8 ).^^,!^")}] , (1) 



where py\ s is the transverse momentum of a "visible object" from a squark/gluino decay, 
which is defined as the sum of visible particle momenta. The p™ lss is the total missing 
transverse momentum. The minimization is taken with respected to the unknown LSP 
momenta p™} ss , p™2 SS under the constraint p™ ss + p™ 2 ss = p™ lss . The transverse mass, m^, 
is defined as 

__,2 /„vis _miss\ f,™Vis\2 , 2 , r> / r-ivis rpmiss _vis _miss\ (r)\ 

™t [Pn > Pt; ) ~ K m i ) +m x + 2 \h Ti h Ti - p Ti ■ p Ti ), (2) 



where E Ti = yp^i + m^- It should be noted that the true LSP mass (m^o) is unlikely to 
be known in advance, so rriT2 is regarded as a function of a test LSP mass (m x ). 

One of the important features is that the is smaller than the parent gluino/squark 
masses if a test LSP mass is set equal to the true value. 

^T2(^ x ?) — max ( m <j? m g)- (3) 

From the upper end point of vtiti (m™| x ), one can obtain the information on the mass 
of the parent particle. Without knowledge of the true LSP mass, m™| x provides a one- 
dimensional constraint between the masses of the squark/gluino and the LSP. 

Recently, it is pointed out that the m™2 X ( m x) function has a kink structure at which 
m x is the true LSP mass unless the squark/gluino decays directly into the LSP through 
a two body decay. An analytic expression of m™| x is derived in Refs. [TT| [T2] . If one 
considers events in which the squark and the gluino are produced in pair with a vanishing 
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total transverse momentum, m m | x (m x ) is given as follows. 



where 



max/ \ \ F< % f° r m Y < m Y° /.\ 

J ^ v XJ \ JF™ ax (m x ) for m x > m x o, v y 



^TK) = -^K is = <L^2 is = <in^ = o,m x ), 

^K) = JK s = ml,mr = <L^ = 0,%). (5) 



Here the function T is given in Ref. [12] and mj is is kinematically bounded as follows, 

m mm < mf < m m L. (6) 



Notice that the events at the end point satisfy mj m = m^ n for m x < m x o while mj m = 
m max f° r m x > m x? - ^ ne kink structure of m™| x (m x ) appears since the functional form 
of m™2 X ( m x) changes at m x = m x o. In Ref. [TT], it is shown that the kink structure 
appears even if the pair-produced squark and gluino have a non-vanishing transverse 
momentum. If one can identify the position of the kink from the LHC experiment, one 
can determine the masses of the squark/gluino and the LSP simultaneously. In Ref. [T2] . 
it is demonstrated that the masses of the squark/gluino and the LSP are determined 
using exclusive decay channel by performing Monte Carlo simulations. In particular, for 
the case that the gluino decay g — > qqxi occurs through the off-shell squark exchange 
diagram, the m™2 X from the gluino pair production has a very sharp kink structure and 
the masses are determined precisely. 



3 Inclusive rriT2 Analysis 

3.1 Hemisphere analysis and inclusive m^2 parameter 

In this section we argue that the kink method discussed in Secf2] can be extended to an 
inclusive analysis. In case of exclusive analyses, one needs to specify a cascade decay 
chain. The branching ratio of the cascade decay chain would depend on the model pa- 
rameters. On the other hand, inclusive distributions are rather insensitive to branching 
ratios. Therefore, if an inclusive quantity can be defined, it may be useful to determine 
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the squark and the gluino masses in the early stage of the LHC experiment. One disad- 
vantage of inclusive approaches may be that all the production and decay modes should 
be taken into account simultaneously, and distribution may be smeared. 

To define an inclusive txiti distribution, we group the final particles into two "visible 
objects". For this purpose, we adopt a hemisphere method in Ref. [HI [H]. For each 
event, two hemispheres are defined and high pt jets, leptons, and photons are assigned 
into one of the hemispheres as follows; 

1. Each hemisphere is defined by an axis p] m (i = 1,2), which is the sum of the momenta 
of high px objects belonging to hemisphere i. We require px > 50 GeV for jets to 
reduce QCD backgrounds. 

2. High pt objects k belonging to hemisphere i satisfy the following conditions: 

d( Pi ,pf)<d(p k ,pf s ), (7) 
where the function d is defined by 

d{ Vk ,pr) = {E i -\p™\cos9^ L: 



Here is the angle between Pi and p k . 

To find axises pj m , we adopted the algorithm discussed in Ref.|13j [T3]. Once p] ls 's are 
determined, one can calculate by using Eq.([T]). 

The inclusive m^ 2 may be compared with MTGEN |15j . The MTGEN variable is a 
minimum of the variable for all possible choices of two subsets of particles a and (3. 
The correct choice of subsets a and (3 leads the heaviest sparticle mass as the end point, 
and the end point should be bounded from above by the mass, if the initial state radiation 
can be ignored. In our algorithm, we assume that the algorithm described above groups 
high pt jets from the same cascade decay with enough probability. This approach is useful 
if the visible hemisphere mass is small compared with the the leading jet energies, which 
is expected especially for the events near the rriT2 end point when the test mass is smaller 
than the LSP mass, see Eqs.(|4]), ([5]). 
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A: MMAM 


B: mSUGRA 




m = 0, R = 20, 
M 3 (GU 1 j = 650 


mo = 1475, mi/2 = 561.2, 
A = 0, tan p = 10 


9 

ul 
Ur 
dR 
X°i 


1491 
1473 
1431 
1415 
487 


1359 
1852 
1831 
1830 
237 



Table 1: Relevant SUSY mass parameters at points A and B. All the mass parameters 
are given in GeV. 

3.2 Model points 

To perform a Monte Carlo analysis, we choose two sample points, A and B. The point 
A corresponds to the MMAM model [TBI T7\ [TBI [19] . In the MMAM model, the mass 
spectrum is parametrized by the modular weights for matter fields n.j and the gravitino 
mass (m 3 / 2 ) and R = m 3 / 2 ((T + T*)/F T ) where T and F T are a modulus field and its 
F-component, respectively. In general, the MMAM model predicts a degenerate SUSY 
spectrum compared with the mSUGRA model. If a = Rj ln(M p i/m 3 /2) is large, the 
SUSY spectrum becomes more degenerate. In this analysis, we choose the point studied 
in Ref . [20] . n.j = 0(1) for squarks and sleptons (Higgs boson), R = 20, tan/3 = 10 and 
the gravitino mass is determined so that M 3 = 650 GeV at the GUT scale. The point 
B corresponds to the mSUGRA with mo = 1475 GeV, min = 561 GeV, A = and 
tan/3 = 10. 

The mass spectrum of SUSY particles is calculated using IS A JET [21J for each sample 
point. In TableUJ, the relevant SUSY masses are listed. At point B, rriq > rrig and 
the gluino undergoes three-body decay through the off-shell squark diagram. The total 
production cross section of SUSY events at the LHC is a = 0.13 pb for both points. 
Squark-gluino coproduction is larger than squark-squark and gluino-gluino productions 
for both points. 

The point B is chosen so that the M e s distribution of one lepton mode is very similar 
to that for point A, where M e g is defined from the sum of the pt of the first four jets and 



6 



a lepton and the missing transverse momentum as follows, 

M eS = j2PT l + PTi + E^ iss . (9) 
i=i 

For the Monte Carlo analysis, we generate 5 x 10 4 SUSY events by HERWIG 6.5 [22] for 
each sample point. To estimate event distributions measured by the LHC detector, we 
use AcerDET [23]. This code provides a simple detector simulation at the LHC. 

In Figfjja), the M e g distribution is shown for one lepton channel. Here we require the 
the following cut. 

1. rijeti^pT > 100 GeV) = n wo > 1 and n- iet {pT > 50 GeV) = n 50 > 4 within \r}\ < 3. 

2. Ef ss > 0.2M efT and Ef ss > 100 GeV and S T > 0.2. 

3. There is one isolated lepton with px > 20 GeV. 

The solid (dashed) histogram is the distribution for point A (B) and the M e g distributions 
roughly agree. 

Although there is not much difference in M e g- distribution defined in Eq.([n]), there 
are more high pt jets on average at point B compared with point A. This is because 
the squark-gluino coproduction is dominant, and a squark decaying into a gluino leads 
additional high pt jets in the events. If one sums all jets with p T >50 GeV, then the 
distribution of point B is significantly higher than that of point A. In FigJT](b), the M e g 
distribution summing up all jets with px > 50 GeV is shown for one lepton channel. We 
will see in the next subsections that inclusive analyses give us a more quantitative 
measure to the difference of the two points. 

3.3 Monte Carlo analysis: Point A (MMAM) 

First, let us consider point A. We require the following cut to select the events. 
1- ""-loo > 2 and n 50 > 4 within \r}\ < 3. 

2. The effective mass of the event must satisfy M eff > 1200 GeV. 

3. At least two jets in each hemisphere. 

4. Ef ss > 0.2M cfT and Ef ss > 100 GeV. 
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Figure 1: (a). M c g distributions for one lepton channel. The solid (dashed) histogram is 
for point A (B). (b). M e g distributions for one lepton channel as in (a) but summing up 
all jets with px > 50 GeV. The solid (dashed) histogram is for point A (B). 

5. There is no isolated lepton with p T > 20 GeV. 

With these cuts, the standard model backgrounds are expected to be reduced significantly, 
so we do not consider the SM background in this simulation. 




400 600 800 1000 1200 1400 

m T g(30) (GeV) 



Figure 2: The two-dimensional distribution in the m^\-R plane for point A. The test LSP 
mass is assumed as m x = 30 GeV. 



To check how well the hemisphere method works, let us consider the following ratio. 

. rn T2 (m x ) - mj%(m x ) 



R( m> 



(10) 



s 



where the parton level raj. 2 is defined so that each visible momentum is the difference of 
a initially produced sparticle and the daughter LSP momentum, and E^ ms is the missing 
energy provided by AcerDET after smearing. If the hemisphere method works well, R ~ 0. 
In FigJ2l the two-dimensional distribution in the ra^\-R plane is shown for m x = 30 GeV. 
The peak of the distribution appears around R ~ 0, but the deviation can be large. The 
reconstructed m^ 2 tends to be smaller than . The main source of the deviation is the 
mis-grouping of the visible objects under the hemisphere method, and also neutrinos and 
jets with pt < 50 GeV which are not included in the hemisphere. 

Let us consider the vtiri distribution for m x < m x o. In Figj3]^a), the m^l distribution 
is shown for m x = 30 GeV. There is an end point at — 1250 GeV. In Figj3](b), 
the reconstructed distribution is shown for m x = 30 GeV. Compared with the m!^\ 
distribution, there is a long tail due to the mis-grouping of the hemisphere method, 
although there is some structure at mri ~ 1250 GeV. 




m T §(30)-30 (GeV) m T2 (30)-30 (GeV) 

Figure 3: (a). The — m x distribution at parton level for m x = 30 GeV. (b). The 
reconstructed uiT2 — w x distribution for m x = 30 GeV. Fitting functions of the end points 
are also shown, see text. 

Let us examine the events around the end point region in detail. FigJH(a) shows 
max(m{ ls , m™) distribution for m T2 (30) > 1000 GeV. As discussed in Section [2j the true 
end point event is realized for mj is = m™ in and the events with large max(m[ ls , m£ s ) are 
considered as fake events. To reduce them, the m^2 distribution for max(m™, m™) < 
400 GeV is plotted in FigJH(b). With the cut on the hemisphere mass, the long tail of 
rriT2 disappears and one can see a rather clear end point at rriT2 ~ 1350 GeV. 
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Figure 4: (a). The distribution of max(m™, m™) for mr 2 (30) > 1000 GeV. (b). The 
mr2(30) distribution for max(m™, rri£ s ) < 400 GeV. 

Next, let us consider the m<r 2 distribution for m x > m x o. In FigJS] (a), the — m x 
distribution is plotted for m x = 900 GeV. There is an end point at ~ 1900 GeV. 
In Figj5] (b), the reconstructed m^2 — fn x distribution is plotted for m x = 900 GeV. The 
distribution has a long tail and one cannot see a clear end point. 




1000 1200 1400 
m T g(900)-900 (GeV) 




1000 1200 1400 
m T2 (900)-900 (GeV) 



Figure 5: (a). The m<T2 — m x distribution at parton level for m x = 900 GeV. (b). The 
reconstructed — m x distribution for m x = 900 GeV. Fitting functions of the end 
points are also shown, see text. 



As discussed in Secj2l the cusp structure of m™^' appears since the functional form 
m™2 X (m x ) changes at m x = m x o. The end point event for m x < m x o is different from 
the one for m x > m x o and these end point events are interchanged at m x 



m v o . To 



confirm this, let us consider how the events near the end point for mr2(30) behaves 
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when m x is large. In Figj6](a), the mT2(900) distribution is plotted for 1200 GeV< 
mr2(30) <1400 GeV. There are two peaks in the distribution. The lower peak is smaller 
than the true end point mT2(900) ~ 1900 GeV. These events are true end point events 
of 7717^(30) while the events around the higher peak are fake events. In FigJHl(b), the 
mT2(900) distribution is plotted using events above the true end point my 2 (30) > 1400 GeV 
We find no peak lower than 1900 GeV as expected, because they are fake events for 
m T2 (30). 




m T2 (900) (GeV) m T2 (900) (GeV) 

Figure 6: (a). The m T2 (900) distribution for 1200 GeV< m T2 (30) < 1400 GeV (left), 
(b). The mT2(900) distribution for the fake events mr 2 (30) > 1400 GeV (right). 

To find the end points of the rriT2 distributions, we also show the fitting of the distri- 
bution in Figsj3] and [51 We fit the reconstructed rrtT2 distribution with a linear function 
which changes the slope at some m x . For comparison, we also show the fitting of rn^\. 
We use a Gaussian smeared fitting function in Ref. [5] for that. xV n -d.f is not good for 
both fits, therefore our fits should be regarded as crude estimates. In addition, the end 
point for m x = 900 GeV depends on the bins used for the fit. Note that the end point 
for m x > m x o is realized for the events with m vls ~ m™^, while the efficiency to assign 
the particles correctly in hemisphere should be low in such case, see FigfTl 

In FigJHl the end points of for various test LSP masses are plotted with solid lines. 
The end points of the m^2 are larger than m^l by 150 — 200 GeV, one can see a kink 
structure around m x ~ 400 GeV, which is close to the true LSP mass, m x o = 487 GeV. 
At the kink, the end point value of is m T2 ~ 1650 GeV. It should be noted that the 
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a) b) 

Figure 7: Kinematical configurations for (a). m vls ~ ui^ n an d 3 )- ^ V1S ~ m max- When 
m vls is large, jets in the hemisphere are less collinear, and the hemisphere analysis likely 
misgroups the particles. 

inclusive mri distribution is dominated by the events from the squark-gluino coproduction 
since the production cross section is larger than those of gluino-gluino and squark-squark 
pair production. In such a situation, the end point of m,T2 distributions is sensitive to 
max(m§, nig). At point A, the gluino is heavier than the squarks and the end point should 
be sensitive to the gluino mass, m g = 1491 GeV. While the end point value is larger than 
the true gluino mass by about 150 GeV, we think the agreement between rriT2 and 
is reasonable given the crudeness of our fit. 

In the mSUGRA, the bino-like LSP mass is about 1/6 of the gluino mass. If we take 
the measured w,T2 end point at the kink as the gluino mass then the LSP mass assuming 
mSUGRA is around 270 GeV. The observed kink is clearly above 270 GeV, therefore we 
can say that the mass spectrum is the MMAM type for this case. 
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Figure 8: The end point of mT2(%) — itl x f° r var i° us test LSP masses. The solid line is 
the m™| x while the dashed line is the parton level my] max . 
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3.4 Monte Carlo analysis: Point B (mSUGRA) 

At point B, the squark is much heavier than the gluino and our interest is to measure the 
squark mass scale quantitatively using the inclusive m,T2 distributions. We require the 
following cuts to select the events. 

1- ^100 > 2 and > 6 within \rj\ < 3. 

2. M eS > 1500 GeV 

3. At least two jets in each hemisphere. 

4. Ef ss > 0.2M cfT and Ef ss > 100 GeV. 

5. There is no isolated lepton with > 20 GeV. 

In Fig|9j the two-dimensional distribution in the m^\-R plane is plotted for m x = 
30 GeV. The peak of the distribution appears R ~ —0.2, and the misreconstruction rate 
is higher than point A. The low reconstruction efficiency may be understood as follows. At 
this point rrig ~ 1850 GeV, m g ~ 1360 GeV and the squark decaying into the gluino gives 
high pt jets as discussed earlier. When this jet is misidentified to the other hemisphere, 
the reconstructed may become much lower than the expected value, because 
the squark is so much heavier than gluino. It can be as low as the order of gluino mass. 
Note that (rrig — rrig) /rrig ~ 0.26, roughly corresponds to the observed shift. Luckily, 
strongly peaks near the end point, and the reconstructed event still makes visible end 
points. 

In Fig JTOT a). the rn^\ distribution is plotted for m x = 30 GeV. There are two peaks 
in the m^l distribution. The higher peak corresponds to the squark while the lower peak 
corresponds to the gluino. At point B, the squark is much heavier than the gluino, so the 
end point is determined by the squark decay, rn^\ ~ 1850 GeV. In Fig JToT b). the rriT2 
distribution is shown for m x = 30 GeV. The witi distribution is smeared but one can still 
see the end point. The end point events are dominated by the events with small m ms . 
There is again interchange of the events near the end point as we increase the test LSP 
mass, and we can see the two peak structure in mT2{ m x > m °xi) distribution for the event 
near the end point of m^^x < m vi) distribution, similar in FiglHl 
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Figure 9: The two-dimensional distribution in m^-ii plane for point B. The test LSP 
mass is assumed as m x = 30 GeV. 




m T g(30) (GeV) m T2 (30) (GeV) 

Figure 10: (a). The m!^\ distribution for m x = 30 GeV. (b). The rriT2 distribution for 
m x = 30 GeV. 

In FigJTTl the end points of m^2 for various test LSP masses are plotted with a solid 
line. The end points are determined as in point A. The end points of the rriT2 is almost 
the same as the ones of m!^\ within errors. One cannot see a clear kink structure around 
the true LSP mass, m x o = 237 GeV. While it is difficult to determine the squark mass 
from the kink method, the inclusive analysis is useful to obtain the information on 
the squark mass. To see whether the end point of m<r 2 correctly describes the squark 
mass for rriq > nig, we show the end point for m x = 30 GeV at the mSUGRA points 
where the gaugino mass is kept the same as that of point B but the universal scalar mass 
mo is varied. In FigJT2| we plot the m<r 2 end point as the function of the squark mass and 
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find very good agreement from 1400 GeV to 1800 GeV. 
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Figure 11: The end point of mT2(%) ~ m x f° r var i°us test LSP masses. The solid line is 
the m™2 X while the dashed line is the parton level my] max . 
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Figure 12: The end point of rriT2 for m x = 30 at several mSUGRA points where the 
gaugino masses is the same as that of point B but is varied. Here the horizontal axis 
is the left-handed up-type squark mass. 
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4 Summary and Conclusion 



In this paper, we have proposed an inclusive analysis at the LHC to obtain informa- 
tion on squark and gluino masses by the hemisphere method. The hemisphere method is 
an algorithm to group collinear and high p T particles and jets, assuming that there are 
two of such groups in a event. The algorithm is to group the cascade decay products into 
two visible objects. To study the distributions, we perform the Monte Carlo simulation 
for two sample SUSY spectra from the MMAM and the mSUGRA models. The cas- 
cade decay products from a squark/gluino are grouped into a visible object with enough 
probability to see the parton level rriT2 end point. However, the end point of the 
distribution is sometime smeared by mis-identification of hemispheres, which obscure the 
end point determination. We have fitted the ttit2 distribution near the end points. The 
end point determination suffers from various systematic uncertainties, such as a choice of 
the fitting function and the fitting region. 

We have examined the events near the end point in detail. For m x < m x o the end point 
events has the minimum hemisphere mass, mj is = m^f n . The fake end point events due to 
the mis-grouping of the hemisphere are reduced if we impose the cut on the hemisphere 
mass, without disturbing the correct end points. For m x > m x o, the end point event is 
realized when the hemisphere mass is at maximum m] ls = m™ x . The end point event 
is interchanged at m x = m x o and the cusp structure of m™| x appears. By checking 
the test mass behavior of the rriT2 variable for the events near the end point, we can 
prove if events near the end point obtained by the fit of iriT2 distribution is correctly 
reconstructed ones or not. We have shown that the true end point event for m x < m x o 
gives the my 2 value smaller than m™% x for m x ^> m x o, while the fake end point event for 
m x < m x o gives larger than x for m x > m x o . From this observation, while there 
are various uncertainties for the end point determination, we conclude that the inclusive 
m T2 distribution is useful to obtain the information of the masses of the gluino/squark 
and the LSP at the LHC experiment. 

For both of the sample points, the main QCD production process of SUSY particles is 
squark-gluino coproduction, and the end point of distribution is sensitive to max(mj, 
trig). For the sample point in the MMAM model, rrig > rriq and the end point should be 
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determined by rrig, while it should be determined by rrig for the sample point in the 
mSUGRA model because rrig > m g . From the Monte Carlo analysis, we have found that 
the end point is indeed determined by max(mj, rriq) for both of the sample points. For 
the MMAM sample point, we have found that there is a cusp-like structure of m™| x (m x ) 
at the true LSP mass and we can determine the gluino and the LSP mass simultaneously. 
For the mSUGRA sample points with rrig > rrig, we find the linearity between rrig and 
m™| x . We have checked that the squark mass is reconstructed up to rrig ~ \Am g when 
m~g ~ 1.4 TeV. 

There have been different approaches in the LHC physics study. One of the direction 
is to study the inclusive quantities such as M eff , E™ 1SS , which do not require reconstruction 
and is useful to grab characters of the events. The other direction is to study quantities, 
which are specific to some processes, such as the end point measurements of the invariant 
mass and the m^2 distribution. They are very powerful to determine the absolute sparticle 
masses. In this paper, we propose an inclusive rriT2, which is inclusive in the sense that 
we do not specify the decay channel. However, it bites the merit of exclusive analyses 
with helps of the new understandings of the stransverse mass function my 2 (m x ) • While 
detailed analyses on systematical uncertainties are still needed, we hope that this quantity 
helps to determine the sparticle masses at the early stage of the LHC experiments. 
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